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ABSTRACT 



Context. The source of emission from Sgr A*, the supermassive black hole at the Galactic Center, is still unknown. Flares and data 
from multiwavelength campaigns provide important clues about the nature of Sgr A* itself. 

Aims. Here we attempt to constrain the physical origin of the broadband emission and the radio flares from Sgr A*. 
Methods. We developed a time-dependent jet model, which for the first time allows one to compare the model predictions with flare 
data from Sgr A*. Taking into account relevant cooling mechanisms, we calculate the frequency-dependent time lags and photosphere 
size expected in the jet model. The predicted lags and sizes are then compared with recent observations. 

Results. Both the observed time lags and size-frequency relationships are reproduced well by the model. The combined timing and 
structural information strongly constrain the speed of the outflow to be mildly relativistic, and the radio flares are likely to be caused 
by a transient increase in the matter channelled into the jets. The model also predicts light curves and structural information at other 
wavelengths which could be tested by observations in the near future. 

Conclusions. We show that a time-dependent relativistic jet model can successfully reproduce: (1) the quiescent broadband spectral 
energy distribution of Sgr A*, (2) the observed 22 and 43 GHz light curve morphologies and time lags, and (3) the frequency-size 
relationship. The results suggest that the observed emission at radio frequencies from Sgr A* is most easily explained by a stratified, 
optically thick, mildly relativistic jet outflow. Frequency-dependent measurements of time-lags and intrinsic source size provide strong 
constraints on the bulk motion of the jet plasma. 
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1. Introduction 

The compact radio source Sagittarius A* (hereafter S gr A*), at 
an estimated distan ce of 8.4+0.6 k pc (Reid et al. 2009) and mass 
of ~ 4 x 10 6 M Q dGhez et al . 2008; Gillessen et al. 2009), is our 
closest supermassive black hole. When compared to other ac- 
tive galactic nuclei (AGN), Sgr A* is extremely under-luminous 
at all wavelengths, suggesting either that the mass accretion 
rate (M) onto the black hole, or the radiative efficiency, is very 
low. Polarization m easurements suggest that M < lO~ 7 M /yr 
(Ma rrone et alj[2007b . Recent multiwavelength campaigns that 
have monitored Sgr A* simultaneously from radio to X-rays 
have found that even though the luminosity of Sgr A* is low 
relative to its AGN counterparts, its emission is quite variable 
at all wavelengths. Flaring activity at both sub-mm and higher 
frequencies has been observed on timescales ranging from min- 
utes to hours, suggesting that the emission at these frequencies 
is produced very close to the supermassive black hole, probably 
within a few tens of gravitational radii (Rg-GM/c 2 ). Although 
flarin g appears to be nearly simultaneous at these frequen- 
cies jMarrone et alj|20fjl lEckart et aLll2008t iDodds-Eden et all 
2009), the flares seem to have an inc reasingly long tim e -dela 
at longer wavelengths. For example, lYusef-Zadeh et al] j200 
found that the 43 GHz flares precede the 22 GHz flares by 20+10 
minutes, and that 350 GHz flares precede those of 43 GHz by 
30 - 60 minutes. 



The physical origin of the observed electromagnetic ra- 
diation from Sgr A* has been under active debate for more 
than a decade. Attempts to model the emission from Sgr A* 
can be broadl y class i fied into two categories: (1 ) accretion in- 
flow models jMeliaHp H iNaravan et all Il998t iLiu & Melial 
l2002h lYuan et al.ll2006l). and (2) relativistic outflow/jet models 
(IFalcke & Biermannlll999t IFalcke & Markofdl2000t lYuan et al.1 
2002). The radio, NIR, and X-ray flare light curves were 
mode led assuming expansion o f a spherical plasm a blob (see 
e.g.. lYusef-Zadeh et all 120081: lEckart et all 12009b . Recently 
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Falcke_etaLll2009) suggested that the radio time-lag data and 
measurements of intrinsic size at various frequencies (e.g., 
iBower et al J [20041: IShen et al.ll2Q05t iDoeleman et ai1l2008l) can 
be used to resolve the degeneracy between inflow and outflow 
models. It h as long been suspec t ed that radio flares such as those 
reported by Yusef-Zade h et al.l ([2008) coul d be "smoking gun 
evidence" of jets in sources such as Sgr A* (Falcke 1999b). 

While jets have not been directly imaged in Sgr A*, sev- 
eral observations sttongly suggest the presence of jets in this 
source. The flat/slightly inverted radio spectral energy disttibu- 
tion (SED) (si gnature of a compact sel f-absorbed jet in the ra- 
dio; see e.g. jBlandford & Konigl|[l979) observed from Sgr A* 
are similar to those of other low-luminosity AGNs. M81 and 
Sgr A* have very similar spectra and polarization properties in 
the cm-radio band , and M81 has been found to have weak jets 
(Nag ar et al J 12005). Rad io observations o f the stellar X-ray bi- 
nary system A0620-00 (IGallo et al.ll2006l) at almost near quies- 
cent luminosity suggest that compact jets are present at bolo- 
metric luminosities as low as ~10~ 7 Lndd- If jet physics and 
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accretion flow scale with the mass of the black hole, then Sgr 
A* (Lt,oi~lCr 9 Lgdd) is also expected to harbor a faint, com- 
pact jet. Steady-state jet models have previously been used to 
model the broadband quiescent spectra o f Sgr A* successfully 
dFalcke & MarkoflfcOOOtlYuan et al.l2002l) ; steady-state models, 
however, cannot predict the flaring activity. 

In this paper we develop a proper treatment of the time- 
dependent cooling processes of the particle distribution, while 
keeping the b asic hydrodynamic outflow model the same as in 
Falcke (1996), allowing a semi-analytical treatment of the prob- 
lem. Once dominant cooling mechanisms are taken into account, 
the time-dependent jet model naturally accounts for all the ob- 
servational constraints, namely, (a) the broadband SED, (b) the 
observed light-curve morphologies and time-lags, and (c) the ob- 
served frequency-size relationship. 



2. Model 

We assume that a fraction of the accretion inflow is chanelled 
into two symmetric, collimated supersonic outflows. The en- 
ergy distribution of the radiating leptons entering the jet is as- 
sumed to be a thermal, relativistic Maxwell-Juttner distribution. 
Most of the kinetic energy of the jet is assumed to be carried 
by cold, non-radiating baryons. The bulk speed at the base of 
the jet is assumed to be the sound speed, and beyond the base 
(or the "nozzle") the jet plasma accelerates longitudinally via 
pressure gradient and expands laterally with the sound speed 
B s ~ 0.4. We adopt the standard one-flui d approximation used i n 
magnetohydrodynamics (MHD; see e.g. iLeismann et al1 l2005). 
and treat the magnetized jet plasma as a single-component fluid 
of adiabatic index 4/3. The longitudinal velocity profile is ob- 
tained by solving the relativistic Euler equation, and th e mag- 
netic field strength estimated via flux conservation (see iFalcka 
1 1 996t iFalcke & Markofll20l)0l for details of the hydrodynamics 
of a pressure driven supersonic jet). 

The observed spectrum from the compact jet is modeled as 
the sum of emission from cylindrical segments along the jet axis. 
Within each segment we consider the following processes influ- 
encing the local evolution of the particle distribution: 

- Synchrotron cooling: Losses due to synchrotron emission are 
given by y Swi =Aca T U bJ 2 /(3m e c 2 ), where cr T is the Thomson 
cross-section and Ub-B 2 /87r is the magnetic energy density. 

- Inverse Compton (IC) cooling: Losses due to IC scatter- 
ings are given by % m =4c(cr m U rad (t)^y 2 /(3m e c 2 ), where, 

following IdeKool et all d 1989b we define (cr KN U rad (t)) = 

4n/c § cr KN (y,v) J Y dv. Here o~ Ktf (v,y) is the Klein-Nishina 
correction to the scattering cross-section and J v the mean in- 
tensity. 

- Adiabatic expansion: Cooling due to adiabatic expansion is 
given by y o ,,=(l/3)y V • v, where V ■ v den otes the divergence 
of the bulk velocity field jLongair II1992I) . 

For numerical discretization we assume that the plasma in 
each segment cools for the time it requires to cross the segment. 
We follow the evolution of such a "parcel" of jet plasma as it 
moves from one segment to another, i.e. in a comoving frame, 
without any particle loss/escape. For simplicity we assume that 
this comoving "parcel" of particles cools and radiates indepen- 
dently, i.e. is not influenced by its neighbors. While this is obvi- 
ously a simplification, given the overall simplicity of the model 



we feel this is justified. The continuity equation for the time evo- 
lution of leptons in this case is given by 

dN(y t) d 

df ' = J- [f (% t)N(y, f)] ; y(y, t) = y Sv „ + y fo ,„ + y id (1) 

Given the low luminosity of Sgr A* (~ 10~ 9 LnddX the com- 
pactness parameter (£=LcrT/(Rnie^y, see Guilbert et al. 1983) 
is very small ({~10~ 6 for R=10 R g ). Therefore pair processes 
are not important here. Also, there is no injection term in ([TJ 
since partic le acceleration seems to be very weak or lacking 
in Sgr A* dMarkoff 1 120051) . Recasting dD as a Fokker-Planck 
equation, we used the f ully implicit Chang-Cooper algorithm 
(IChang & Cooperill970t [Chiaberge & Ghisellini 1999) to solve 
it. 

Once the time-evolved particle distribution solution is ob- 
tained, we then compute the emitted SE D due to angle averaged 
synch rotron emission using eq.(10) of ICrusius & Schlickeiserl 
(Il986h for relativi stic electrons with Loren tz factor >2, and 
using eq. (8) of iKatarzvri ski et alj d2006l) for electrons in 
the cyclo-synchrotron regime with Lorentz factor <2. IC is 
computed incorporating the Klein-Nishina correction for scat- 
tering beyond Thomson li mit following the prescription of 
iBlumenthal & Gould! dl970t) . The seed photons for IC are the 
photons produced locally by synchrotron emission (synchrotron 
self-Compton; SSC). In the absence of detailed spatial informa- 
tion of the hot star s that constitue the central star cluster close 
to Sgr A* (see e.g. iGenzel et al.ll2000T) . we estimated emission 
from these stars by assuming 100 stars of temperature 10,000 K 
each, situated 0.5 parsecs from Sgr A*. The photon density of 
locally produced synchrotron+SSC photons in the jet is at least 
3 orders of magnitude higher than that due to total photon den- 
sity from these stars, suggesting that external inverse Compton 
losses due to photons from the nearby cluster is not important. 

Since the bulk motion of the emitting jet plasma is relativis- 
tic we first compute the Doppler factor T> = [T(l - Bcosff)Y l , 
where 9 is the angle between velocity and line of sight, cB is the 
bulk speed, and F = 1 / - B 2 is the bulk Lorentz factor. Then 
the observed SED is calculated by applying the appropriate spe- 
cial relati vistic transformations of the emitted frequency and flux 
(see e.g. lLind"& Blandf ordl 19851) . and taking into account time- 
lags due to the spatial extent of the flow. 

For any given frequency, when the flux from each jet seg- 
ment is plotted as a function of distance from the black hole, it 
roughly resembles a bell-shaped curve (in the self-absorbed part 
of the jet). The full width at half maximum (FWHM) is an indi- 
cator of the size of the photosphere of the jet at this frequency. 
We used this procedure to calculate the frequency-size relation- 
ship. 

In our model the radio flares are caused by a perturbation 
in M at the jet-base, leading to a transient density enhancement. 
This density enhancement during a flare of duration tj is quan- 
tified by the parameter f„ so that the number density at the base 
(«o) as a function of time is given by 



«o(0 = 



(1 +/„) n (0 < t < t d ) 
«o otherwise. 



(2) 



As the overdensity propagates outward along the jet, the ad- 
ditional pressure in the overdense region would cool and expand 
into the neighboring segments with lower pressure and density. 
While a full computation of this diffusion (and resultant cool- 
ing) would require MHD simulations beyond the scope of the 
present paper, the effect can still be captured by assuming an ad- 
ditional velocity dependent cooling term in the continuity equa- 
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tion for the overdense region. Assuming that the overdense re- 
gion expands longitudinally into the neighboring segments with 
some fraction of the sound speed, we write this additional cool- 
ing term as a function of distance from the black hole (z) as 
yjy = f o p s c/z- Thus f f3 s c can be interpreted as the longitudinal 
speed with which the overdense region expands. 



3. Modeling the quiescent broadband emission and 
the radio flares from Sgr A*on July 17, 2006 

We test the model against the simultan eous observations of Sg r 
A* in 22 and 43 GHz on 2006 July 17 dYusef-Zadeh et al.ll2008h . 
which has the best simultaneous coverage to date of a flare in 
both frequencies. Indirect evid ence suggests that Sgr A* is seen 
under large inclination angle (Markoff et al. 2007; Falck eetafl 
2009). For simplicity we therefore assume that the jets are per- 
pendicular to our line-of-sight, even though the model can han- 
dle any inclination angle. The free parameters of the model are: 
temperature (T e ) and density (no) of the thermal leptons at the 
base, location of the sonic point or the nozzle (ho), radius of the 
nozzle (ro) and its magnetic field (Bo), The flare is parametrized 
by its start time (to), duration (td), density enhancement fraction 
(f n ) and expansion speed of the overdense region (f , as a frac- 
tion of the sound speed). 

The good agreement of the data with model light curves is 
shown in Fig. Q] The model parameters are given in Table Q] 
Note that we only model the flare which peaks around 6.5 hours 
UT in 43 GHz. The increase in flux in both bands during the end 
of the obs ervations may be th e beg inning of another flare as sug- 
gested by Yus ef-Zadeh et alj d2008l) . but we do not model this 
due to lack of full coverage. In Fig. [2] we show model predicted 
light curves at 0.7, 1.3, 2, 3.6, 6 and 13 cm. From Fig. Q]it is clear 
that the light curves are asymmetric and even a phenomenologi- 
cal model would require (a) rise time-scale, (b) decay time-scale, 
and (c) amplitude. In our model also three parmeters describe the 
flare completely, viz. tj (which can be roughly associated with 
the rise rime), f (associated with decay time), and f n (associ- 
ated with amplitude). It is thus important that a single set of flare 
parameters fits two rather disparate light curves at two frequen- 
cies. Moreover, the decrease of amplitude is consistent with the 
decrease in average variability with frequency. A compilation 
of the rms flux variab ility of Sgr A* averaged o ver many flares 
bvlFalckeetalJ (|2009|) from VLA and GBI data dFalckdll999at 
iHerrnstein et alj 20041) shows that the rms variability decreases 
with increasing wavelength. The relative amplitude of this flare 
at both 43 and 22 GHz is about 0.6 times the long-term rms vari- 
ability, therefore this flare was not extraordinarily bright or faint, 
and may be considered as a fairly typical flare from Sgr A*. 

The frequency-size relationship predicted by the model is 
displayed in Fig. [3] showing reasonable agreement with the data 
at low frequencies where the jet becomes optically thick. The 
discrepancy at shorter wavelengths where the spectrum is opti- 
cally thin is expected given that our model ignores general rel- 
ativistic effects, and simplifies the physics of jet formation and 
launching, both of which must play an important role very close 
to the black hole. For this model, the rate at which matter is fed 
to both jets is M = 2nom p c/3 s y s nrQ = 2.4 x 10~ 9 M e yr" 1 , which 
is about 2 orders of magnitude lower than the upper limit of 10~ 7 
M Q yr _1 suggested by linear polarization observations. 



4. Discussion 

In this paper we have shown that a time-dependent relativistic 
jet model can explain most of the observational features seen in 
Sgr A*. The main conclusions of this work are: 

- The model presented here can describe the quiescent broad 
band SED of Sgr A*, as well as its long-term radio variabil- 
ity. 

- Assuming that the radio flares are caused by a transient den- 
sity enhancement at the base of the jet, which then prop- 
agates outward, the model can also fit the 22 and 43 GHz 
light curve morphology o f the flares seen on 2006 July 17 by 
Yuse f-Zadeh et alj (2008). We predict light curve morphol- 
ogy and delays at longer wavelengths, which can be used to 
test the model with future observations. 

- The model predicted frequency-size relationship also 
matches quite well with the observed radio data. 

The radio emission in our model originates in the self- 
absorbed, optically thick part of the outflow, where adiabatic 
losses dominate over radiative cooling. Therefore variability in 
the radio light curves largely traces the expansion of the jet 
plasma. Assuming that the flares have the same frequency-size 
relationship as the quiescent emission, the time delay measure- 
ments combined with frequency-size measurements strongly 
constrain the bulk speed of the jet plasma in this model and 
rule out subrelativistic expansion speeds. An alternate model 
used by lYusef-Zadeh et al] d2008l) to fit the radio light curves 
assumes adiabatic, subre lativistic expansion o f a spherical blob 
of plasma, i.e. the classic Ivan der Laanl d 19661) model. Our mod- 
eling, at the very least, shows that the van der Laan model is 
not unique and in fact the jet model, in contrast to the van der 
Laan model, not only correctly predicts the light curves and the 
overall variability amplitudes, but also fits the frequency-size 
relation without violating any other observational constraints. 
It must be kept in mind that the jet model does not attempt 
to explain the NIR or X-ray flares. The NIR/X-ray flares trace 
particle (re)energization and cooling very close to the black 
hole, while only margin ally affecting the optically thick radio 
flux (Mar koff et alj|200ll) . and no tight correlation between ra- 
dio and X-ray flares have been found so far. A full general 
relativistic MHD model including radiative coo ling (see e.g. 
iFragile & Meieil 120091: iMoscibrodzka et al.|[2QQ9h would be re- 
quired to model the jet launching region. However the simple 
model presented here captures the important physics and shows 
that the jet model explains the radio flaring properties naturally, 
due to simple adiabatic expansion of overdensities in the outflow 
(likely linked to variations in the accretion rate). Future coordi- 
nated multiwavelength campaigns, especially measurements of 
time-lags, sizes and positional offset at other wavelengths, will 
enable a better understanding of the velocity profile of the jet 
from Sgr A*. 
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Table 1. Model parameters for the flare of 2006 July 17. The 
mass (Mb h) and distance were taken from their latest estimates 
(see e.g. iGillessen et al] 2009: Re id et al.l 12009). Inclination of 
the jet axis to the line of sight was fixed to 90 degrees. The base 
of the jet is parametrized by radius of the nozzle (rn), location of 
the sonic point (ho), magnetic field strength (Bo), number den- 
sity (no) and temperature (T e ) of the thermal leptons. The flare is 
characterized by its start time (to), duration (to), density enhance- 
ment fraction (f n ), expansion speed of the overdense region (f a , 
as a fraction of the sound speed). Parameters marked with an 
asterisk (*) were not varied during fitting. 
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Fig. 1. Comparison of model with data for the flare on 2006 
July 1 7 at 43 and 22 GHz. T he 43 GHz data with errorbars 
(from lYusef-Z adeh et al. 2008) are shown in red, and 22 GHz 
data+model in blue. As noted in |3] we model the flare which 
created a peak in the 43 GHz light curve near 6.5 hours UT. 
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Fig. 2. Model predicted light curves of Sgr A* at key radio fre- 
quencies using model parameters shown in Table Q] Also shown 
as error bars between t=200-300s are the rms flux variability of 
Sgr A * averaged over many flares as compiled by iFalcke et al.l 
(2009). The gray bands have a width of 0.6 times the rms vari- 
ability of Sgr A*, which is roughly matched to the specific flare 
under consideration. 
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Fig. 3. Comparison of model predicted frequency-size re- 
lation ship (solid line) with obser v ations . The data are 
from [Bower et all d2004l) : IShen et al.l d2005l) : iDoeleman etail 
(200 8]), which wer e compiled and corrected for scattering by 
IFalcke et all (120091) . 
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